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Abstract 

We present the measurement of the lifetime of the SI 6 we have performed using 
three independent data samples from two different decay modes. Using a T,~ beam of 
340 GeV/c we have obtained clean signals for the SI 6 decaying into S - A' - 7r + 7r + and 
Sl _ 7T + 7r _ 7r + , avoiding topological cuts normally used in charm analysis. The short but 
measurable lifetime of the SI 6 is demonstrated by a clear enhancement of the signals 
at short but finite decay lengths. Using a continuous maximum likelihood method 
we determined the lifetime to be rno=55T 3 ( s ^ a C723( s y s C f s - This makes the SI 6 the 
shortest living weakly decaying particle observed so far. The short value of the lifetime 
confirms the predicted pattern of the charmed baryon lifetimes and demonstrates that 
the strong interaction plays a vital role in the lifetimes of charmed hadrons. 
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1 Introduction 


The experimentally observed lifetimes of charmed baryons differ by as much as a factor of 
three between r(E+) = 350t][o f" s [l]-[5] and r(E°) = 98t 23 fs [6, 7, 5]. This large variance 
is understood to be the consequence of the strong interaction between the quarks involved 
and of W-exchange between the charm quark and a d-quark [8, 9]. 

If both decay modes are possible, the direct decay of the c-quark by emission of a W + 
(c —► sud) and W-exchange with a d-quark of the baryon (cd —► su), the lifetime is short, 
as is the case for the . If only direct decay is possible, the lifetime is long, as is the 
case for E+ . Further differences among the lifetimes follow from the details of the baryon 
wavefunctions. 

The decay is of particular interest. It proceeds dominantly via the direct decay. In 
the final state, however, three identical quarks are present (0°—► sssud ). If the two strange 
quarks of the are only spectators, the lifetime will be comparable to the E+ lifetime. 
If, however, the three strange quarks overlap in phase space in the final state, the decay 
will be accelerated because there are three ways to reach the same final state, and the 

lifetime will be very short. In addition the Cabibbo suppressed W-exchange with an 
s-quark (cs su) may contribute [10]. 

The has been observed so far in three different experiments: in the hyperon beam ex¬ 
periment WA62 [11], in the e + e _ experiment ARGUS [12] and the photon beam experiment 
E687 [13, 14]. A preliminary estimate of its lifetime has been reported by the collaboration 
E687 at the HQ94 conference, r (fl°) = 50-100 fs [15], superseded by a result presented 
at LISHEP 95, r (D°)= 86±^ fs [16]. In the following, we present a measurement of the 

lifetime from the hyperon beam experiment WA89 at CERN, using data recorded in 
1993. Preliminary results were already presented at Moriond 95 [17]. 

2 Beam and apparatus 

The hyperon beam is derived from a 450 GeV/c proton beam with an intensity of 
4-10 10 p/spill (2.5 seconds every 14.5 seconds) which interacts in a 40 cm long beryllium 
target. Negative particles are extracted in the forward direction by a magnetic beam chan¬ 
nel of 12 m length defined by a tungsten collimator. Downstream of the collimator, beam 
particles are detected in a set of 9 silicon microstrip detectors of 50 /im pitch and two 3 mm 
thick scintillators used in the trigger. The beam has an intensity of 6-10 5 particles/spill and 
is composed of tt~ and E _ with a ratio of about 2.3:1 at the experiment target. The mean 
momentum of the E _ is 340 GeV/c. An 80 cm long transition radiation detector (TRD) 
serves as a fast tagger on E _ [18]. The signal of the TRD is used in the 2nd level trigger to 
reject about 95% of all high energy beam 7r - , keeping 85% of the E _ . 

The apparatus is designed as a forward spectrometer and schematically shown in Fig 1. 
It consists of a target area, a decay zone and a magnetic spectrometer followed by detectors 
for particle identification and calorimetry. 

The beam impinges on a longitudinally segmented target composed of one copper slab 
(4 mm thick) and 3 carbon plates, made of pressed diamond powder (p ~ 3.3 g/cm 3 ) with 
a thickness of 2 mm each. The targets are placed in a closed box flushed with helium and 
are spaced by 2 cm. A vertex detector made from 29 planes of silicon detectors is placed 
downstream of the target box. The first 24 planes have a readout pitch of 25 pm and active 
areas of 5x5 cm 2 and 6.3x6.3 cm 2 [19]. The last 5 planes have an active area of 5x5 cm 2 
and a readout pitch of 50 pm. Ten detectors measure the horizontal x-coordinate, 10 the 
vertical y one. Six detectors measure the coordinate (y+45°) and 3 the coordinate (y—45°). 
Two 3 mm thick scintillators mounted behind the 12th silicon plane are used to trigger on 
interactions, thus the first part of the vertex detector also serves as silicon targets of 3.6 
mm total thickness. 

The decay zone downstream of the target area has a length of about 10 m. It contains 
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drift chambers (40 planes) which are interleaved with 1 mm pitch MWPCs with an active 
area of about 12x12 cm 2 including the beam region (20 planes). Their purpose is to detect 
track segments from hyperon and K° decays before they enter the magnetic spectrometer. 
MWPCs with 1 mm pitch (12 planes) placed 2 m downstream of the target are used to 
link tracks reconstructed separately in the two parts of the detector. The particle momenta 
and charges are determined in the Omega spectrometer, a large superconducting magnet 
(f B ■ dl ~ 7.2 Tm) equipped with MWPCs (45 planes) measuring three different co¬ 
ordinates (horizontal and ±10.15°). Two sets of large size drift chambers (12 planes) are 
used to measure the tracks of particles leaving the spectrometer. They are mounted back 
to back with MWPCs (MY1 and MY2) used for pattern recognition and triggering. 

A large Ring Imaging Cherenkov (RICH) detector serves for the identification of particles 
with momenta above 12 GeV/c [20]. It provides 7 t/K separation up to 100 GeV/c. Two 
hodoscopes, HI and H2, are mounted on either side of the RICH and serve as trigger 
detectors. Downstream of H2 we use a lead glass calorimeter with an active area of 3 m 2 
[21] followed by a hadron calorimeter made from lead/scintillating fibres constructed in the 
spaghetti calorimeter technique [22]. 

The main trigger consists of three levels. At the first level, a beam trigger and a signal 
from the two interaction counters corresponding to two or more particles are required. At 
the second level, signals from > 3 charged particles are required in HI and MY1 and hit 
correlations between HI, H2 and MY1 are used to enrich the sample with events having 
at least two charged particles above 30 GeV/c. E _ interactions are selected at the third 
level using the TRD (see above). This rather loose interaction trigger and the geometrical 
acceptance of the spectrometer limit the acceptance of the experiment to charmed hadrons 
produced with Xj > 0.15 — 0.2, depending on the final state. 

Track reconstruction is performed separately in each detector area. Track segments are 
then linked towards the upstream direction, from the spectrometer to the decay zone and 
then to the vertex detectors. Track segments in the vertex detector linked to track segments 
in the spectrometer are refitted taking into account multiple Coulomb scattering. Beam 
tracks are reconstructed in the beam detectors and are required to match the geometry of 
the magnetic beam channel to reject particles from hyperon decays or from interactions in 
the walls of the collimator. 

3 Selection and Reconstruction of 0° Decays 

The data sample comprises 1.5-10 8 interactions on tape which correspond to about 2 • 10 10 
incoming E _ and 3.8-10 8 interactions. The decay channels used in this analysis are 

fl°^fl - 7r + 7r - 7r + (1) 

£l 0 c ^E~ K~ir + ir + (2) 

The selection of candidates for charm baryon decays is based on particle identification, 
reconstruction of the interaction point and isolation of the secondary vertex. 

A—► p7r — are identified through their decay geometry and kinematics, A K~ and 

E - ^- A7T - decays are reconstructed from negative particles forming a vertex with a A and 
giving an invariant mass near the nominal mass of fl - or E _ . Only combinations are kept in 
which the momentum sum points to a track in the vertex detector (assumed to be the 0“ or 
E _ ) which has no matching track in downstream detectors. In order to reduce background 
in the 0“ sample we exclude track combinations with a A7T - mass consistent with a E _ . 
For positively RICH-identihed K~ with momentum greater 25 GeV/c the window around 
the E _ mass is ±5 MeV/c 2 , whereas below that value the window is ± 10 MeV/c 2 and we 
require that in addition the tt~ hypothesis is rejected by the RICH if the K~ candidate 
falls into its geometrical acceptance. The invariant mass plots for the hyperons used in this 
analysis are shown in Fig. 2. 

The reconstruction of the charm decay vertex follows a ’candidate-driven’ approach. 
First we select track combinations for which the particle types or at least their charges 
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correspond to charm decays. A common vertex then is reconstructed for these tracks (’sec¬ 
ondary vertex’). From the remaining tracks we determine the position of the ’main vertex’, 
iteratively omitting tracks from the combination until the quality criterium (x 2 /d/ <5) is 
met [23]. At least three tracks are required for a main vertex. 

In the following we describe the selection of the charm candidates in more detail. The 
criteria used are generally the same for all samples allowing small variations of their values 
depending on the final state. Reconstructed secondary vertices are retained if they fulfill 
requirements on the y 2 of the vertex (x 2 /d/ <5-6). We also require that the calculated 
position error for the main vertex is below 300-500 [im and below 1 mm for the charm 
decay vertex. Tracks with momenta below 1 — 2 GeV/c are not considered for secondary 
vertices. 

In order to obtain clean mass distributions, charm samples are generally selected using 
also other criteria, not used for the final lifetime analysis. These are based on the topologies 
of these two vertices (e.g. vertex isolation and individual track impacts) as well as require¬ 
ments on a spatial separation r usually expressed in terms of its significance r/oy. This 
significance is calculated from the position errors of the two vertices which are determined 
from the error matrices of the tracks used. Fig. 3 shows the invariant mass spectra obtained 
for 3 different data samples using a significance cut of about 3 a. These comprise the decay 
channels fi _ 7r + 7r _ 7r + (henceforth referred to as sample A) and K~tt + 'it + . The latter 
channel has been split into two independent samples, one using only K _ well identified 
in the RICH (sample C), the other one using the complementary sample (with K~ below 
the Cherenkov threshold or outside the RICH acceptance) coming from interactions in the 
carbon target (sample B). The reason for the exclusion of interactions in the other targets 
from the last sample will be explained below. Although decays have been reconstructed 
in more decay channels [24, 25] here we only present those final states used for the lifetime 
measurement, selected by statistical significance and life time resolution. Table 1 gives a 
summary of the signals from Fig. 3. 


channel 

mass [MeV/c 2 ] 

width [MeV/c 2 ] 

S/B 

A) fi _ 7r + 7r _ 7r + 

2707.5±3.9 

8.3T2.9 

23/18 

B) S - K~ 7r + 7r + (Carbon) 

2720.0±3.4 

8.8T3.4 

62/162 

C) S - A - 7r + 7r + (high K-momenta) 

2719.8±4.7 

8.9T3.0 

27/50 


Table 1: Properties of the signals from Fig. 3 obtained from a likelihood fit to the 
data. The errors quoted are statistical only. Systematic errors in particular on the mass 
determinations are currently under study. The observed differences in mass for the three 
measurements have, however, no influence on the determination of the lifetime. S/B is 
given for a mass windows of ±15 MeV/c 2 . 

Topological selections such as the 3 a cut mentioned above usually cause distortions of 
the decay distributions which are not easy to represent by Monte Carlo simulations. This 
is particularly true for very short lived particles. It therefore is the aim of this analysis to 
determine the lifetime of the without such cuts. 

In the following we present the evidence for a short but measurable lifetime of the 
using all three data samples. Fig. 4a shows the distribution of the measured decay 
times 


t — (r sec ^prim) j C 

for events off the signal region (see Fig. 3) from sample A. No minimal separation between 
main and secondary vertex has been required. The distribution reflects the measurement 
accuracy obtained and can be fitted empirically using two Gaussians with widths of 31 
and 69.5 fs with relative weights of 1:2. They are well reproduced in the simulation. In 
Fig. 4b we show the same distributions for events from the signal band which shows a 
shift towards positive values indicating a short but measurable lifetime for the Fig. 5a 
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shows the same distribution for sample B. Again, the distribution can be described by two 
Gaussians of width 48 and 125 fs and relative weights of about 4:1. Fig. 5b gives the mass 
distribution for these combinations using a separation cut of 3 a for positive and negative 
vertex separations. We clearly see that a signal only exists for positive (physical) lifetimes 
and that the background is similiar in both cases at the level of 5-10%. This again is clear 
evidence for a finite lifetime. Finally we show in Fig. 6 the invariant mass distributions for 
the sample C which contains K~tt + 'it + combinations with a kaon positively identified in 
the RICH. The distributions are plotted for different minimal decay lengths required of 0, 
50, 100 and 150 fs. The short but measurable lifetime of the is clearly seen. For this data 
sample the resolution on the decay times is again described by two Gaussians (<Ji= 45 fs, 
(T 2 = 105 fs) with relative weights of 4:1. 

The origin of the different values for the time resolution of different data samples lies in 
the different underlying momentum spectra of the candidates. Sample B exhibiting the 
longest tails shows the lowest average momentum. 

Samples with low energy kaons from interactions in copper have not been considered. 
Due to the large contribution of a component with a =150 fs no clear signal could be 
obtained without additional, topological selection criteria such as vertex isolation. This 
wide distribution can be explained by the large amount of multiple scattering to which the 
low momentum secondaries are subjected since the copper target is the most upstream one. 


4 Determination of the Lifetime 


For the determination of the lifetime, a sample is chosen from a mass window around 
the peak centre (’signal band’). A second sample from mass windows away from the peak 
([2500,2650] MeV/c 2 and [2750,2900] MeV/c 2 ) is used to study the background (’side band’). 
The mass window chosen for the signal band is ±15 MeV/c 2 corresponding to about ±2<r 
obtained from the fits to the mass spectra. The lifetime is determined from a maximum like¬ 
lihood fit to the decay time (/) distribution of the signal band using a continuous likelihood 
method. 

The likelihood function for each event consists of an exponential with lifetime r convo¬ 
luted with a resolution function /„(/) and a function f b {t) which describes the background: 

1 r°° 

£r(<) = rijy-J o f.(t-0)e-'' T dO + (1 ~v)fb(t) 

with the normalisation factor 

pOO poo 

N t = / f s (t-0)e~ s/T d0dt 

Jo 


depending on the cut on the minimum accepted decay time. In this analysis we assumed 
that the time distribution of the background events reflects the time resolution and therefore 
we set f s (t)=f b (t). 

The resolution functions 





D -t 2 /(2<±) 


1 — ot j 2 / / *-> 2 \ 

_ e ~ t /C 2 ^) 

o 2 V2vt 


consist of two Gaussians for which the parameters are obtained from a fit to the side band 
distribution (see above). The background events in the signal band are assumed to have 
the same distribution as in the side bands. 

The parameter tj, giving the fraction of the signal in the mass window considered, is 
estimated from a fit to the observed mass distributions. Its value is varied in the lifetime fit 
under the constraint of a Poisson distribution which is determined by the number of events 
in signal and side bands. 
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The lit results are shown in Fig. 7. They are displayed as a function of the lowest 
lifetime t min still used in the fit. The errors are the statistical errors of the fits. The results 
for the different samples agree very well over the full range of t min . The fit results are also 
shown in Table 2. The average values for the lifetime are obtained from a simultaneous 
likelihood fit to all three samples. 

Table 2: Lifetimes measured in different samples as a function of the lowest decay time t min 
used in the fit. 


sample 

^ min [^] 

T [fs] 

V 

Combinations in 





signal band 

+ 

1 

+ 

1 

all 

qo +15 
oz _ 12 

.27 ±.08 

112 


0 

oq +15 
zo _ 14 

.41 ±.10 

68 


25 

26 ti® 

.49 ±.11 

43 


50 

2Q + 25 
-22 

.45 ±.14 

26 


75 

50 til 

.36 ±.18 

11 


100 

63 t 9 ° 

.39 ±.24 

4 

+ 

+ 

1 

1 

[I] 

all 

29 

fixed at .08 

1158 

from Carbon 

0 

20 t\l 

.18 ±.04 

594 


25 

36 ±i 5 5 

.21 ±.04 

374 


50 

57 ±i| 

.25 ±.05 

207 


75 

63 t\l 

.30 ±.06 

128 


100 

74 til 

.34 ±.07 

64 

+ 

+ 

1 

1 

[I] 

all 

58 + 2 ( 

.14 ±.05 

259 

K ident. in RICH 

0 

qo +30 
OO _28 

.25 ±.07 

151 


25 

59 til 

.29 ±.08 

93 


50 

62 tH 

.40 ±.09 

56 


75 

70 tit 

.45 ±.11 

34 


100 

oq +50 
OO _ 28 

.48 ±.14 

20 

all channels 

all 

28 tl 




0 

25 tl 




25 

37 




50 

55 t\l 




75 

64 t] 5 2 




100 

77 ±i 9 5 




Studies of systematic errors 

In the following we discuss various sources of systematic errors. The strongest systematic 
influence on the value of the lifetime comes from the variation of the lower decay time cut (see 
Fig. 7). We have investigated this effect with simulated events passing the same selection 
criteria. A similar systematic behaviour was observed, however of smaller amplitude. We 
also observe a systematic downward shift of the reconstructed main vertex which, however, 
shows little influence on r for values of t min > 50 fs. These shifts can be compensated by 
using Gaussians not centred at the origin. However, we do not correct our data for this 
effect but attribute a systematic error of Ar=28 fs, allowing for the lowest and highest 
lifetime value observed to be included at their la level. 

Another source of uncertainty comes from the parametrisation of the resolution function. 
For each decay it is chosen to be the sum of two Gaussians fitted to the time distributions 
in the sidebands (see above). Since this function serves two purposes, as indicated by 
distinguishing between /„(/) and /+), we have investigated them separately. If the wider 
Gaussian is omitted completely in the resolution function /„(/) for the signal events, but 
retained for /+), the lifetime increases by 7-9 fs for t min in the range of 0 to 100 fs. The 
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Gaussians do not fully describe the background time distribution at t > 200 fs for the 
fi _ 7r + 7r _ 7r + mode and at t > 300 fs for the K~tt + 'it + mode. A bias from this should 
show up in a variation of the lifetime if an upper limit t max is applied to the fitted time 
spectra. For our chosen value t min = 50 fs, the lifetime varies by ±10 fs if t max is varied 
from ±oo down to below 200 fs. We therefore attribute a total systematic uncertainty of 
1 jo fs t° account for possible errors connected with the resolution function or large timings. 
Secondly we have investigated the influence of the description of the background by 
Variations in the fraction of the long living component by ±10-20% cause variations in the 
fitted lifetime, in particular for small values of t min . At large values ( t min > 50 fs) this 
effect becomes small (±5 fs). The sample most affected is sample B which exhibits the 
largest tails in the resolution function. The variations are of the order of the statistical 
error for t min < 25 fs. However, the influence on the averaged lifetime is small (±4 fs). The 
sample fi - 7r + 7r - 7r + , on the other hand, shows only little dependence (±5 fs for variations of 
about 20% on the fraction of the long living component), as the resolution function exhibits 
shorter tails. The same holds for sample C where only the data point with t min = 0 shows 
some significant variation. The total systematic error attributed to uncertainties in fi{t) is 
±5 fs. 

Since the resolution for secondary vertices depends on the momentum of the daughter 
particles we have checked possible biases stemming from the momentum spectrum. We have 
subdivided our samples rejecting about 30% of the events from low and high energy 
candidates, respectively, and repeated our analysis. We have found a small systematic 
dependence on the momentum of the order of ±5 fs. 

In addition, we have investigated the sensitivity of the fits to the choice of the mass 
window chosen for the signal. Only a small variation (±5 fs) is observed varying the mass 
window between ±10 and ±20 MeV/c 2 . 

From the present understanding of our detector, from our reconstruction algorithms and 
analysis cuts we do not observe a significant influence on the lifetime: biases of the resulting 
reconstruction efficiencies are small thus we omit the application of a correction function to 
our measured decay time distributions. Finally we investigated a possible dependence on 
the choice of the background windows chosen as well as on the y 2 cuts used for the secondary 
vertices. No significant influence was observed. In particular the timing distributions for 
the background windows below and above the mass are alike. 

Since the influence of systematic effects becomes smaller for larger lifetime cuts ( t min > 
50 fs) we choose the data point at t m in= 50 fs as central value since it has the best statistical 
significance. This choice leads to an asymmetric systematic error evaluation due to the 
dependence of r on t min . The total systematic error is estimated to be fs. 

5 Conclusion 

We have measured the lifetime of the obtaining the final result: 

r(fi°)= 55tjj(sfa±)tjgjsr/sf.) fs. 

This value makes the the shortest living weakly decaying particle observed so far. The 
short lifetime has to be attributed to the properties of the final state containing three 
strange quarks and shows that the values of charmed hadron lifetimes depend largely on 
the structure of the baryon. 
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Combinations/2 MeV/c' 


WA89 

Hyperon beam 
1993 layout 



DC,MWPC MWPC DC Hod 1 Hod 2 

Target area Decay area Omega spectrometer RICH Calorimeters 


Beam Si TRD Beam Si Targets Vertex Si 



Scin Scin 12 x 25micron 6.5x6.5cm 
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5 x 50micron 5.0x5.0cm 
Target area Angles: 0,90,-45,+45 


Figure 1: WA89 experimental set-up in 1993 
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Figure 2: Examples for mass distributions of reconstructed strange particles. 
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Figure 3: Mass distributions for three different final states: a) fi - 7r + 7r - 7r + , b) 
"“A'“7r + 7r+ from carbon, c) S _ A’~7r + 7r + from all targets with positively RICH identified 
kauri. The shaded region denotes the signal band. 
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Figure 4: ft 7r + 7r 7r + sample : vertex separation (ps) for a) background events, b) for signal 
baud. 


13 



























































































































vertex separation (ps) 


o 

> 

CD 


.2 

c5 

c 

£ 


o 

u 



Figure 5: K~Tr + Tr + sample in carbon: a) vertex separation for background events (ps) 

b) " - A' - 7r + 7T + mass spectrum after requiring 3<r separation of the two vertices (negative 
separations: filled histogram, positive separations: open histogram). 
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Figure 6: invariant mass distributions for the H FT 7r + 7r + sample with positively identified 
kaons using different lifetime cuts for the charm candidates. 
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Figure 7: Results of the likelihood fits (see text) for the different samples and their depen¬ 
dence on the lower decay time cut t min used in the fit. The shaded area represents the size 
of the full systematic error. Sample A) denotes l) - 7r + 7r - 7r + , sample B) S“A’~7r + 7r + from 
the carbon target and sample C) K~Tr + Tr + with identified K~. 
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